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ABSTRACT 


A discrete  vortex  method  has  been  formulated  for  calculating 
the  hydrodynamic  loads  on  SES  sidewalls  in  an  inviscid  separated 
flow  due  to  vorticity.  A rigid,  flat,  free  surface  is  assumed. 
The  vortex  system  consists  of  a bound  vortex  lattice  attached  on 
the  lifting  surface,  a single  free  tip  vortex  core  from  the  side 
edge,  and  free  vortex  lines  emanating  from  the  trailing  edge. 

The  presence  of  these  free  vortices  accounts  for  the  nonlinear 
effects.  Single  and  double  sidewalls  of  different  area  aspect 
ratios  have  been  considered  for  various  geometrical  arrangements 
and  flow  conditions.  It  was  found  that  with  the  present  single 
tip  vortex  model,  the  accuracy  of  the  results  improved  with  in- 
creasing aspect  ratio.  Side  forces  predicted  for  small  AR 
cases  were  higher  than  given  by  other  methods.  The  effect  of 
wall  interference  was  significant  for  double  sidewalls  with  in- 
creases in  length-to-beam  ratio,  or  area  ratio,  or  yaw  angle. 


ADMINISTRATIVE  INFORMATION 

This  work  was  authorized  by  the  Naval  Sea  Systems  Command  (SEA  035) 
and  funded  through  the  General  Hydromechanics  Research  Program  under 
Task  SR0230101,  Work  Unit  1-1606-040.  The  study  was  conducted  by  Codes 
1606  and  1615  and  monitored  by  Code  1506,  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center  (DTNSRDC). 

INTRODUCTION 

Recent  improvements  in  surface  effect  ships  (SES)  have  been  aimed  at 
higher  performance  operations  of  SES's  with  high  length-to-beam  ratios. 

The  sidewalls  of  such  crafts  are  subject  to  severe  side  forces  and 
moments  which  sometimes  become  detrimental  to  maneuverability.  Acting 
as  low-aspect-ratio  lifting  bodies,  the  hydrodynamic  loads  on  those 

; 

sidewalls  are  greatly  influenced  by  the  presence  of  vortices  created  by 
flow  separation  from  the  leading  and  side  edges  of  the  surface. 

The  problem  is  essentially  nonlinear  because  of  the  free  vortex 
influence.  Even  though  flow  separation  exists,  the  flow  field  can  still 
be  analyzed  by  using  an  inviscid  model  when  separation  occurs  only 
along  the  sharp  edges.  Limited  theoretical  work  has  already  been  performed 


in  this  area.  Whitehead  developed  a simplified  model  by  assuming  a 

rigid  free  surface  on  the  upper  ends  of  the  sidewalls  and  neglecting  the 

2 

effects  of  cushion  pressure  and  seals.  Lamar  used  a generalized 

suction  analogy  to  predict  force  contributions  due  to  the  vortices  shed 

from  the  sidewalls.  In  conjunction  with  the  side-edge  suction  analogy, 

3 

Woodward  et  al.  obtained  the  required  pressure  distribution  from  a 

separate  three-dimensional,  lifting  surface  theory  of  attached  potential 

flow.  They  found  reasonable  results  for  rectangular  parallel  sidewalls 

with  aspect  ratios  of  0.1  and  0.2,  but  did  not  determine  the  location  of 

the  shedding  vortices  because  of  the  limitation  of  the  side-edge  suction 

analogy.  This  same  limitation  precluded  calculating  the  influence  of  one 

sidewall  when  a vortex  is  shed  from  an  adjacent  sidewall.  The  use  of 

linear  potential  flow  theory  also  raises  some  questions  for  cases  with 

aspect  ratios  of  0.1  or  less. 

4 

Wilson  also  employed  the  generalized  suction  analogy  but  he  cal- 
culated the  required  pressure  distributions  from  slender  wing  theory 
rather  than  linear  potential  flow  theory.  For  most  cases,  his  results 
on  net  side  force  and  moment  compared  fairly  well  with  those  obtained  by 
Whitehead. ^ However  the  two  approaches  gave  different  trends  for  force 
and  moment  coefficients  on  each  individual  wall  as  affected  by  interference 
due  to  spacing. 

The  present  work  considers  the  problem  of  two  parallel  sidewalls  at 
a constant  yaw  angle  and  utilizes  a discrete  vortex  method  to  simulate 


Whitehead,  R.E.,  "A  Theoretical  Prediction  Technique  for  Lateral 
Hydrodynamic  Loads  on  Sidewalls  of  Surface  Effect  Ships,"  DTNSRDC 
Report  4766  (Dec  1975).  A complete  list  of  references  is  given  on 
page  55  • 

^Lamar,  J.E. , "Extension  of  Leading-Edge  Suction  Analogy  to  Wings 
with  Separated  Flow  Around  the  Side  Edges  at  Subsonic  Speeds,"  NASA 
TR  428  (Oct  1974). 

3 

Woodward,  R.A.  et  al.,  "A  Computer  Program  for  Three-Dimensional 
Lifting  Bodies  in  Subsonic  Inviscid  Flow,"  United  States  Army  Air 
Mobility  R&D  Laboratory  Report  TR-74-18  (Apr  1974). 
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Wilson,  M.B.,  "Interference  Effects  on  Lateral  Forces  and  Moments 
on  High  L/B  SES  Arrangements,"  Paper  76-859  presented  at  the  AIAA/SNAME 
Advanced  Marine  Vehicles  Conference,  Arlington,  Virginia  (Sep  1976). 
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the  lifting  surface  of  the  sidewalls  and  the  shedding  of  vortices  from  the 
side  edge  of  the  sidewall  surface.  The  nonlinear  effect  of  flow  separation 
from  sharp  edges  is  accounted  for  by  considering  the  free  tip  and  wake  vor- 
tex cores  away  from  the  surface.  For  simplicity  the  tip  vortices  are 
represented  by  a single  free-vortex  core  which  is  merged  by  side  edge 
vortices  a short  distance  away  from  the  wall.  An  iterative  procedure  is 
employed  to  obtain  the  direction  of  the  core  as  well  as  the  wake  vortices 
in  the  solution. 


VORTEX  FLOW  MODELS 

SINGLE-WALL  MODEL 

First,  we  consider  the  problem  of  an  SES  with  a flat  wall  of  finite 
aspect  ratio  traveling  in  calm  water.  The  upper  end  of  the  wall  is  aligned 
with  the  surface  of  the  water.  With  the  assumption  of  a rigid  free  sur- 
face, the  flow  can  be  simulated  by  a bound  vortex  sheet  and  free  vortex 
sheets  emanating  from  the  side  and  trailing  edges  of  the  wall.  A schematic 
of  the  flow  is  shown  in  Figure  1.  The  yaw  angle  a is  measured  with  respect 
to  the  x-z  plane.  The  z = 0 plane  is  assumed  to  be  flat  and  rigid,  and  its 
influence  is  represented  by  an  image  in  the  (-z)  region. 

The  vortex  system  of  the  subject  boundary  value  problem  consists  of 
a bound  vortex  sheet  attached  on  the  lifting  surface,  a free  edge  vortex 
sheet  rolling  up  from  the  side  edge,  and  a free  wake  vortex  sheet  emanating 
from  the  trailing  edge;  see  Figure  1.  There  is  a pressure  difference 
across  a bound  vortex  sheet  but  not  across  a free  vortex  sheet.  However, 
the  normal  component  of  velocity  vanishes  on  either  a bound  or  free 
vortex  sheet. 

The  present  single-wall  arrangement  is  similar  to  that  of  Kandil^ 
except  that  the  tip  vortex  sheets  have  been  simplified  to  a single  vortex 
core  and  that  no  vortex  segment  is  placed  along  the  side  edge  in  order  to 
satisfy  the  Kutta  condition  there. 


Kandil,  O.A.,  "Prediction  of  the  Steady  Aerodynamic  Loads  on  Lifting 
Surfaces  Having  Sharp-Edge  Separation,"  Ph.D.  Thesis,  Virginia  Polytechnic 
Institute  and  State  University,  Blacksburg,  Virginia  (Dec  1974). 
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DOUBLE-SIDEWALL  MODEL 

In  practical  SES's,  there  are  always  two  sidewalls  placed  side  by  side 
and  separated  by  a finite  distance.  To  simplify  the  problem,  it  is  assumed 
that  the  effect  of  the  cushion  pressure  and  seals  can  be  neglected.  The 
problem  is  then  reduced  to  two  distinct  single  walls  with  mutual  inter- 
ference. The  schematic  of  the  vortex  system  for  the  two-sidewall  model 
is  shown  in  Figure  2.  As  for  the  single-wall  case,  it  is  assumed  that 
the  free  surface  is  rigid  and  flat  and  that  its  influence  is  accounted 
for  by  the  vortex  images  in  the  (-z)  region. 

Mutual  interference  results  from  the  presence  of  a vortex  field  in 
both  sidewalls.  The  velocity  at  a control  point  has  to  be  calculated  by 
considering  the  contributions  of  the  vortex  elements  not  only  from  a single 
sidewall  but  also  from  those  of  the  adjacent  sidewall  a finite  distance 
away.  As  a general  rule,  however,  only  those  vortex  elements  that  the  con- 
trol point  can  "see"  make  contributions.  Because  of  the  presence  of  the 
other  wall,  the  shape  of  the  vortex  sheets  of  the  upstream  sidewall  is 
also  confined  to  limited  space  between  the  solid  boundaries;  however, 
this  assumption  is  important  only  in  the  case  of  very  narrow  spacing  and 
high  yaw  angles. 


BASIC  EQUATIONS 

INDUCED  VELOCITY 

For  an  inviscid  and  incompressible  fluid  flow,  the  velocity  generated 
by  a vortex  segment  is  determined  by  using  the  Biot-Savart  law;  see 
Figure  3.  For  a control  point  o^  in  the  presence  of  vortex  segment  qs 
with  circulation  y,  the  induced  velocity  is  given  by 

V » -~^r  (cos  9.  + cos  6 ) e ... 

4nh  1 2 v (I) 


where 


cos  6 ^ 


(2) 


Here  e is  the  unit  vector  of  the  induced  velocity  and  e is  the 
v qs 

unit  vector  parallel  to  the  vortex  segment  qs. 


The  derivation  of  Equation  (1)  can  be  found  in  standard  text.  As  shown 
in  Figure  4,  a control  point  is  used  to  apply  the  boundary  condition  only 
when  necessary  to  determine  the  influence  coefficients,  as  discussed  later 
in  the  section  on  numerical  procedure. 

INDUCED  FORCE 

The  force  induced  by  a bound  vortex  segment  in  a vortex  field  can  be 
calculated  by  using  the  Kutta-Joukowski  theory.  With  reference  to  Figure  5, 
the  induced  force  on  vortex  segment  ab  in  the  presence  of  vortex  segment  qs 
is  given  by 

F = PL  r v x e (4) 

Here  L is  the  length,  e is  the  unit  vector,  and  F is  the  circulation 
of  segment  ab.  V is  given  in  Equation  (1).  It  should  be  noted  that 
F differs  from  Y of  segment  qs. 

Equations  (1)  and  (4)  are  the  basic  equations  for  computing  the 
necessary  flow  characteristics  of  an  inviscid  and  incompressible  flow. 

They  originate  from  the  nonlinear  partial  differential  equations  govern- 
ing a rotational  flow,  with  the  aid  of  the  theorem  of  Stokes.**  Flow 
separation  caused  by  vortex  rather  than  viscosity  may  take  place.  Such 
a situation  may  exist  in  a flow  with  angle  of  yaw  up  to  26  degrees,  as 
observed  experimentally. 

At  the  surface  of  bound  vortex  sheets,  the  normal  component  of  velocity 
vanishes  and  there  is  a pressure  difference  across  the  sheet;  the  boundary 
conditions  are: 


**Prandtl,  L.  and  O.E.  Tietjens,  "Fundamentals  of  Hydro-  and  Aerodynamics," 
translated  by  L.  Rosenhead,  Dover  Publications,  New  York  (1957). 
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Discrete  Vortex  System  for  Determining  Circulation 


n . V - 0 


AC  0 
P 


At  the  trailing  and  side  edges,  however. 


in  order  to  satisfy  the  Kutta  condition.  The  normal  component  of  velocity 
also  vanisher  at  the  surface  of  a free  vortex  sheet  but  there  is  no  pressure 
difference  across  a sheet  and  so  the  boundary  conditions  are: 

n . V » 0 


C ■=  0 
P 


DISCRETE  VORTEX  METHOD 

The  circulation  values  in  Equations  (l)  and  (A)  are  generally  unknown 
and  must  be  determined  numerically  in  order  to  calculate  th#  velocity  and 
forces.  For  illustrative  purposes,  we  consider  the  single-wall  model  first. 

In  the  numerical  procedure,  a bound  vortex  sheet  is  represented  by  a 
lattice  and  a free  vortex  sheet  by  a series  of  nonintersecting  lines.  For 
simplicity,  those  nonintersecting  free  edge  vortex  lines  that  emanate  at 
the  end  of  each  row  of  the  vortex  lattice  are  merged  into  a single  vortex 
core  at  each  side  of  the  wall;  see  Figure  4.  Both  the  free  wake  vortex 
lines  and  the  free  tip  vortex  cores  are  composed  of  a series  of  short 
vortex  segments,  except  for  the  last  segment  which  is  semi-infinite  and 
extends  downstream  parallel  to  the  free  stream.  These  segments  are 
generally  separated  from  the  surface  and  their  final  location  is  determined 
by  the  solution  with  the  aid  of  an  iteration  procedure.  The  force-free 
condition  on  the  free  vortex  sheet  is  approximately  satisfied  by  aligning 
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these  segments  with  local  streamlines.  The  presence  of  these  free  vortex 
segments  accounts  for  the  nonlinear  relationship  between  yaw  angle  and 
force  coefficients. 

The  lattice  of  the  bound  vortex  divides  the  lifting  surface  into 
panels.  Each  panel  has  a spanwise  vortex  segment  along  its  forward  edge 
and  a chordwise  vortex  segment  along  its  inboard  edge.  The  center  of  the 
panel  is  the  control  point  at  which  the  boundary  condition  of  the  vanish- 
ing normal  velocity  component  n . V = 0 applies.  No  vortex  segment 
is  placed  along  the  trailing  and  side  edges  in  order  to  satisfy  the  Kutta 
condition. 

DETERMINATION  OF  CIRCULATION  FIELD 

When  a cartesian  coordinate  system  is  used,  then  (see  Figure  A)  the 
velocity  at  the  control  point  0^  will  be  the  vector  sum  of  contributions 
from  all  vortex  segments  in  the  field. 


where  N , N , and  N are  the  number  of  vortex  segments  in  the  bound, 
b t w ’ 

tip,  and  wake  vortex  sheets,  respectively.  Subscript  t denotes  spanwise 
segments  and  t,  chordwise  segments.  Terms  with  subscript  i are  contri- 
butions from  the  image.  Other  variables  are: 


(cos 


+ cos  0 ) 

K 


(9) 


‘“‘l  + W 


(10) 


\ ■ *y~  <“*  e«.K  + ~8  v 


(id 


H 9»‘k  + cos  9‘V 


(12) 


g " 7—^ ( COS  0 . + COS  0 - ) 

'M  4,heH  elM  a2« 


(13) 


8ei  ' 4ih~ ' <cos  esll„  + co8  *ol2M> 
M ei^  M M 


(14) 


g - 7~r — (cos  0 . + cos  0 „ ) 

w 4irh  wl„  w2  „ 

N w N N 


(15) 


H ’ <c°’  9-“„ + cos  W 

N 


(16) 


The  variables  0 and  h are  strictly  geometrical  properties  which  can  be 
evaluated  for  a given  configuration  and  vortex  system  arrangement.  The  de- 
tailed evaluation  for  these  quantities  is  given  in  the  appendix. 

However,  the  values  of  circulation  y are  still  unknown  at  this  stage. 
They  will  be  determined  as  follows. 

With  the  aid  of  the  image  principle,  we  write 


Yti  “ Yt 
K K 


(17) 
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From  conservation  of  vortex  at  any  spatial 

v = y 

l 

+ Yt 

K K-N 
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z 

as  illustrated  by  the  following  sketch 
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At  the  edge  of  the  sidewall,  the  circulation  of  the  tip  vortex  line 
at  a given  row  equals  that  of  the  vortex  segment  adjoining  the  tip  vortex 
core,  i.e.. 


y - y.  when  K « M N 
eM  CK 


as  illustrated  in  the  following  sketch 
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At  the  trailing  edge,  the  circulation  of  the  wake  vortex  line  at  a 
given  column  equals  that  of  the  vortex  segment  adjoining  the  wake  vortex 
line  and  thus 


Y - Y when  K > N (N 
\ lK  Z X 


as  illustrated  by  the  following  sketch 


TRAILING f 
EDGE  1 


The  unit  normal  vector  at  control  point  0^  (indexed  by  k for  receiver) 


*n  " 'l  X Vl'l  X ‘2I 


where 

Cj  - (x;J-x1)I  + (y3-y1)j  + (z3~21)k 


C2  " *x4_x2)1  + (y4-y2)J  + (z4_z2)k 


Here  i,  j,  k are  the  base  vectors  in  a cartesian  coordinate  system. 

The  normal  component  of  velocity  at  point  0c  is  thus  obtained  by  sub- 
stituting Equations  (8)  and  (24)  into  (5).  After  simplifying,  we  get 


Gkl  ytt  + Gk2  yt2  + Ck3  yt3  + * * * GkNb  + Rk  “ ° (27) 


where  G^  is  a geometrical  quantity  and  is  the  normal  component  of 
the  free  stream  velocity.  Similarly,  this  equatior  holds  for  other  points 
in  the  bound  vortex  lattice  from  k ■ 1 to  N,  : 

D 
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This  equation  gives  an  N,  number  of  linear  algebraic  equations  for  N. 

b b 

number  of  unknowns  y and  thus  indicates  the  strength  of  the  circulation. 

CK 

Again,  the  influence  coefficients  G are  geometrical  properties.  Several 

KK. 

numerical  procedures  are  available  for  solving  Equation  (28).  The  particu- 
lar one  employed  here  is  the  Gaussian  elimination  scheme  without  pivoting. 
(A  Gauss-Seidel  iteration  scheme  was  also  used  for  some  checkout  cases.) 

The  solution  gives  the  circulation  values  for  y • The  values  for  y , 

K \ 

y ,and  y can  then  be  found  with  the  aid  of  Equations  (21)  and  (23). 
eM  WN 


DETERMINATION  OF  FREE  VORTEX  LINES 

The  geometrical  properties  represented  by  the  influence  coefficients 
depend  strictly  on  the  position  of  the  vortex  segments.  For  free  vortex 
lines  representing  the  free  edge  and  wake  vortex  sheets,  their  positions 
coincide  with  the  local  streamlines  which  are  determined  by  the  solution. 
This  is  done  through  an  iteration  procedure. 

The  positions  of  free  vortex  lines  are  initially  specified  at  the  end 
of  each  row  of  the  bound  vortex  lattice  for  the  free  tip  vortex  lines  and 
at  the  end  of  each  column  for  the  free  wake  vortex  lines.  Since  these 
lines  are  free  of  force,  they  will  eventually  align  with  the  local  stream- 
lines, as  mentioned  previously.  During  the  course  of  solution,  therefore, 
the  positions  of  those  free  vortex  segments  were  updated  by  modifying 
their  coordinates  with  increments  in  the  direction  of  local  velocity. 

The  modified  coordinates  were  then  used  to  determine  a new  set  of 
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geometrically  oriented  influence  coefficients;  subsequently  the  circu- 
lation values  were  obtained  and,  finally,  a new  velocity  field.  The 
coordinates  of  those  free  vortex  segments  were  then  corrected  again  with 
slopes  generated  by  new  velocities.  This  process  was  repeated  until 
further  change  in  coordinates  were  not  in  excess  of  a specified  tolerance. 

CALCULATION  OF  FORCES  AND  MOMENTS 

The  velocity  term  in  Equation  (4)  must  be  determined  before  the  force 
can  be  calculated.  The  procedure  is  similar  to  that  discussed  in  the 
previous  section  except  that  the  control  point  has  to  lie  on  the  vortex 
segment  rather  than  at  the  center  of  a panel,  as  required  by  the  Kutta- 
Joukowski  theory.  Therefore,  with  reference  to  Figure  5,  the  velocity 
induced  at  control  point  0 will  be  the  sum  of  contributions  of  all 
vortex  segments  with  known  circulation 


\ f‘8'K  \ + SilK  YtlK>  VK  + S \ + 8»K  V 


■2X  \ + s 


M-l 


(29) 
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=>wi  Ywi  ^ evt 
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where 


' «hT  (c”  *'ti,  * cos  ’t2> 

N t,.  IV.  K 


(30) 


8ti„  “ W 7 (cos  6^il  + cos  ®ti2  > 
t K.  K. 


8v  ‘ (cos  eiiK + cos  9)iV 

K 


*«„  • ■i^rr  <c°s  **ii„ + coa  **i2.) 

K.  *'1-.  K.  K. 


8e  * (COS  ®e  1 + COS  ®e2  5 

M 4 ^ M M 


g'  = ■ -T7 (COS  0'  + COS  0'.o  ) 

elM  4*hei„  ellM  el2M 


8w  " 4ffh'~  (cos  9^I  + cos  6w2  } 

N wM  N N 

N 


8wi  " (cOS  0wil  + C0S  0wi2  ) 

N wi„  N N 
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Again,  0 • • • and  h • • • 


are  geometrical  quantities  which 


can  be  evaluated  as  before,  but  everything  is  based  on  control  points  0{ 

and  0.  rather  than  0 . 

* c 

The  force  generated  at  the  midpoint  of  the  spanwise  segment  ac  is 
therefore  calculated  by  using  Equation  (4) 


F - pL  r V x e 
0t  t t 0t  vt 
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Similarly,  Che  force  generated  at  the  midpoint  of  the  chordvlse  segment  ab 


F - pL,  T V x e . 

0^  l l 0^  vl 

and  therefore  the  resultant  total  force  on  a particular  panel  K Is 


F = (F  + F ) 

K 0.  ;K 

t * 

It  can  also  be  written  in  a general  form 


F-fi  + fj+fK  (4 

K *K  V ZK 

The  side  force  and  drag  coefficients  are  calculated  by  summing  the 
force  components  respectively  normal  to  and  parallel  to  the  free  stream 


u 

(fyK  cos  ° ■ fx„  sin  a)/(l  p v-2  bc) 


u 

S c<”  “ + £y„  ,ln  *>/ (1  p v-2  bc) 


The  moment  coefficient  with  respect  to  the  leading  edge  is 


\E  -ZWd p v-2b'2) 

L*E*  K-l 


and  that  with  respect  to  the  quarter  chord  is 
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■MHH 


(45) 


(CL  cos  a + Cp  sin  a)  (C.P. 


0.25) 


where  C.P.  is  the  center  of  pressure,  i.e.. 


C.P. 


cos  a + CD  sin  a) 


NUMERICAL  RESULTS  AND  DISCUSSIONS 


(46) 


Numerical  results  have  been  obtained  for  thin  rectangular  single  and 
double  sidewalls  with  virtual  area  aspect  ratios  of  0.05  to  2.0  (or  real 
area  aspect  ratios  of  0.025  to  1.0)  at  yaw  angles  up  to  12  degrees.  For 
all  the  cases,  the  real  plane  contained  40  panels  in  the  bound  vortex 
lattice,  a single  tip  vortex  core  and  four  wake  vortex  lines.  The  tip  vor- 
tex core  was  divided  into  15  segments  and  each  wake  vortex  line,  five  seg- 
ments. The  selection  of  these  numbers  was  a compromise  between  the 
accuracy  and  computer  time.  The  use  of  other  numbers  was  also  investi- 
gated and  typical  results  are  shown  in  Figure  6 for  the  case  of  AR  » 2.0 
and  a « 12  degrees.  The  program  takes  approximately  1 to  1.5  minutes  in 
the  CDC  6600  computer  for  a typical  single-wall  case  and  3 to  4 minutes 
for  a double-wall  case.  Figure  6 also  shows  the  effect  of  initial  location 
of  the  tip  vortex  core  on  its  shape  and  the  side  force  coefficient. 


SINGLE-SIDEWALL  CASES 

Figure  7 shows  the  shape  of  the  tip  vortex  core  and  wake  vortex 
lines  for  a single  wall  of  AR  = 2.0  and  a - 12  degrees.  The  final  lo- 
cations of  these  free  vortex  lines  converged  after  four  iterations.  The 
deflections  shown  were  in  the  x-y  plane  only.  In  the  x-z  plane  which 
is  not  shown,  the  tip  vortex  core  tended  to  sweep  slightly  inward  until 
it  passed  35  percent  of  the  chord  and  then  outward  all  the  way  downstream. 
The  wake  vortex  lines,  on  the  other  hand,  swept  outward  consistently  and 
the  outermost  line  bent  the  most.  The  resulting  pressure  distribution 
is  presented  in  Figure  8.  Note  that  the  leading  edge  suction  was  highest 
close  to  the  center  of  symmetry.  Note  also  that  there  was  a local  positive 


OZ  = HV 


Figure  6 - Effects  of  Number  of  Panels  of  Bound  Vortex  Sheet 
and  Initial  Location  of  Tip  Vortex  Core 


rtex  Core  and  Wake  Vortex  Lines  for  a Single  Wall 


pressure  region  in  the  downstream  area  near  the  center  of  symmetry.  The 
trend  was  consistent  with  that  calculated  by  the  three-dimensional  po- 
tential flow  program.*  The  resultant  normal  force  coefficient  was  in 

fairly  good  agreement  with  experimental  data^  ***  and  with  theoretical 

2 

results  obtained  by  using  the  generalized  suction  analogy;  see  Figure  9a. 
However,  the  present  results  were  fairly  sensitive  to  the  location  of  the 
tip  vortex  core  near  the  leading  edge;  i.e.,  the  farther  away  from  the 
leading  edge,  the  higher  the  side  force. 

The  above  example  with  high  AR  (Figure  9a)  was  considered  for  the 
purpose  of  validating  the  present  numerical  procedures.  More  realistically, 
however,  AR  values  for  practical  SES  sidewalls  are  always  low.  Our  next 
illustration  is  therefore  concerned  with  very  low  aspect  ratio  values  of 
AR  = 0.05.  The  results  for  this  case  at  a * 2,  4,  and  6 degrees  are  de- 
picted in  Figure  9b  along  with  experimental  data  of  Kaplan**  and  theo- 

2 

retical  results  of  Lamar  who  used  the  generalized  suction  analogy. 
Generally,  the  side  forces  predicted  by  both  theories  were  slightly 
higher  than  the  experimental  values.  The  present  calculated  moment  co- 
efficients are  in  good  agreement  with  the  Kaplan  experimental  data. 

The  reason  why  the  present  approach  predicts  higher  forces  in  low  AR 
cases  is  attributed  to  the  limit  of  the  present  single  tip  vortex  core 

^Parkinson,  G.V.  et  al.,  "Observations  on  Low  Aspect  Ratio  Wings  at 
High  Incidence,"  Canadian  Aero  & Space  J. , Vol.  13,  No.  3,  pp.  111-116 
(Mar  1967). 

g 

Winter,  H. , "Flow  Phenomena  on  Plates  and  Airfoils  of  Short  Span," 

NACA  TM-798  (1936). 

9 

Scholz,  N. , "Kraft-und  Druckverteilungsmessungen  an  Tragflachen 
kleiner  Struckung,"  Forsch.  Geb.  Ing.-Wes.,  Vol.  16,  No.  3,  pp.  85-91 
(1949/50). 

***Schoch,  D.L.,  "An  Investigation  of  the  Flow  Characteristics  about 
a Low  Aspect  Ratio,  Sharp  Leading-Edge  Rectangular  Wing,"  Report  574 
(Contract  Nonr  1958(14)),  D.  Aeronautical  Eng.,  Princeton  Univ.  (Nov 
1961);  available  from  DDC  as  AD  270  111. 

**Kaplan,  P. , "An  Experimental  Study  of  SES  Craft  Lateral  Hydro- 

dynamic  Forces  and  Moments,"  Oceanics,  Inc.  Report  73-97  (May  1973). 
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model.  The  inaccuracies  involved  in  using  a single  vortex  core  to  repre- 
sent a sheet  of  vorticity  rolling  up  at  the  sidewall  edge  are  greatest 
for  the  Induced  velocities  in  the  near  field,  i.e.,  near  the  edge  of  the 
wall.  Since  a rectangular  wall  of  low  AR  value  has  an  "edge  region"  that 
is  relatively  much  larger  than  those  of  high  AR  ones,  the  resultant  side 
force  is  therefore  higher  for  the  former. 

DOUBLE-SIDEWALL  CASES 

Several  double-sidewall  cases  were  considered  to  evaluate  the  effect 
of  interference  due  to  various  arrangements.  Area  aspect  ratios  were 
varied  from  0.1  to  0.5  and  yaw  angles  from  2 to  12  degrees.  The  spacing 
between  two  walls  was  arranged  0.1  to  0.3  chord  lengths  apart,  equivalent 
to  length-over-beam  ratios  of  3.33  to  10. 

Figure  10  indicates  the  influence  of  area  aspect  ratio  on  side  force 

at  a = 6 degrees  for  different  length-over-beam  ratios.  The  interference — 

as  indicated  by  the  difference  in  values  for  each  sidewall — became 

more  significant  as  AR  increased.  This  was  also  true  for  narrower  spacing 

between  the  two  walls.  Figure  10b  includes  results  obtained  by  Wilson 

4 

from  slender  wing  theory  with  the  aid  of  generalized  suction  analogy. 

The  comparison  is  rather  inconclusive  except  that  both  indicate  more  side 
force  was  generated  on  the  upstream  sidewall  than  on  the  downstream  one. 

The  values  obtained  by  Wilson  were  based  on  a thickness  of  6 percent;  no 
thickness  was  assumed  in  the  present  approach,  however. 

Figure  11  shows  the  effect  of  length-to-beam  ratio  on  the  side  force 
at  a - 6 degrees  for  an  area  aspect  ratio  AR  = 0.1.  It  is  a crossplot 
of  Figures  10a  and  10b  and  so  the  trend  is  consistent  with  those  figures. 
Note  that  as  the  L/B  value  vanished,  the  side  forces  on  both  walls 
approached  the  single-wall  value.  The  fact  that  the  present  scheme  pre- 
dicts slightly  higher  side  force  than  given  by  slender  wing  theory  is 
attributed  to  the  constraint  of  the  single  tip  vortex  core,  as  already 
discussed. 


lJ 
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Figure  10  - Effect  of  Area  Aspect  Ratio  on  Side  Force 


Figure  10a  - For  L/B  = 5 and  a = 6 Degrees 


10  = uv 


Even  with  interference,  the  linear  relationship  between  side  force 
and  yaw  angle  held  for  small  aspect  ratio  cases  (AR  = 0.1)  with  moderate 
length- to-beam  ratios  (L/B  = 5),  as  observed  in  Figure  12a.  The  same  was 
true  for  the  moment.  For  larger  AR  values  (0.2  to  0.5),  however,  the 
linear  relationship  broke  down,  as  indicated  in  Figures  12b  and  12c  and 
interference  became  more  severe  as  yaw  angle  increased.  The  linearity 
noted  at  the  lowest  aspect  ratio  is  a fortuitous  result  of  the  particu- 
lar aspect  ratio,  length- to-beam  ratio,  and  the  vortex  model.  It  does 
not  indicate  any  linear  asymptotic  character  of  the  theory  for  vanishing 
aspect  ratio. 

The  above  results  were  generated  within  the  capability  of  the  present 
single  tip  vortex  core  model.  The  model  gave  correct  information  on  the 
force  and  moment  characteristics  at  various  flow  conditions  and  geometric 
orientations.  Because  of  the  unavailability  of  adequate  experimental  data 
for  double-sidewall  cases,  it  is  impossible  to  assess  the  interference 
effect  quantitatively  by  the  present  method. 

CONCLUSIONS  AND  RECOMMENDATIONS 

A discrete  vortex  method  has  been  formulated  for  calculating  the 
forces  and  moments  on  SES  sidewalls  in  an  inviscid  separated  flow  due  to 
vorticity.  Single  and  double  SES  sidewalls  with  different  area  aspect 
ratios  were  considered  for  various  geometrical  arrangements  and  flow 
conditions.  It  was  found  that: 

1.  The  present  method  predicts  higher  forces  in  low  aspect  ratio  cases 
than  do  other  methods.  This  is  attributed  to  the  use  of  a single  tip 
vortex  model  to  represent  a sheet  of  vorticity.  The  dominance  of  side- 
wall  edge  effects  reduce  accuracy  at  low  aspect  ratios.  However, 
accuracy  improves  with  increasing  aspect  ratio. 

2.  In  the  case  of  double  sidewalls,  the  effect  of  wall  inter- 
ference becomes  significant  with  increases  in  length-to-beam  ratio,  or 
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Figure  12  - Side  Force  and  Yaw  Moment  Coefficients  at  Various  Yaw  Angles 


Figure  12a  - For  AR  = 0. 1 and  L/B 


Figure  12  (Continued) 


Figure  12c  ~ For  AR  — 0.5  and  L/B 


area  aspect  ratio,  or  yaw  angle.  The  magnitude  of  force  and  moment  for 
either  sidewall  approaches  single-wall  values  as  the  length-to-beam  ratio 
diminishes. 

In  order  to  extend  the  method  to  its  full  capacity,  it  is  recommended 
that  it  be  refined  to  include  multitip  vortex  lines  and  allow  consideration 
of  wall  thickness.  For  verification  purposes,  experimental  investigations 
should  be  conducted  on  double  sidewalls  in  calm  water  with  zero  cushion 
pressure. 


APPENDIX 


CALCULATION  OF  INDUCED  VELOCITY  AND  FORCES  ASSOCIATED 
WITH  DISCRETE  VORTEX  METHOD 

For  a thin  SES  sidewall  traveling  at  a yaw  angle,  a,  and  a pitch  angle, 
8,  the  flow  can  be  simulated  by  a bound  vortex  sheet  on  the  surface  of  the 
sidewall  and  free  vortices  emanating  from  the  trailing  edge  and  from  the 
edge  tip.  The  schematic  of  the  flow  was  shown  In  Figure  1:  a Is  measured 

with  respect  to  the  x-z  plane  and  6 with  respect  to  the  x-y  plane.  The  z-0 
plane  is  assumed  to  be  flat  and  rigid  and  Its  Influence  is  represented  by 
an  Image  in  (-z) . 

A segment  in  a vortex  field  serves  simultaneously  as  sender  and 
receiver.  A sender  induces  flow  in  the  surrounding  region  and  a receiver 
generates  velocity  as  a result  of  induction  from  a sender.  In  a discrete 
vortex  system,  a sender  or  a receiver  may  consist  of  more  than  one  straight 
vortex  segment.  For  convenience,  three  types  of  receivers  and  three  types 
of  senders  are  classified  below;  they  are  similar  to  those  developed  by 
Kandil,^  except  the  change  in  z direction. 


TYPES  OF  RECEIVERS 
Type  1 Receivers 


0 

c 


(n,m) 


(n+l,m) 


This  type  of  receiver  is  found  in  the  bound  vortex  sheet  for  calcu- 
lating circulation  influence  coefficients  with  center  at  control  point 
0c,  whose  coordinates  are  given  by 
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+ yiH-l,m  + yn+l,nH-l^ 

(A-2) 

+ ZrH-l,m  + Zn+l,nrt-l^^ 

(A— 3 ) 

Segment  K in  Figure  4 is  a typical  example  of  a receiver  of  Type  1. 
Velocity  is  computed  at  control  point  0 


Type  2 Receivers 


(n,m+l) 


(n,m) 
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The  velocity  is  calculated  at  midpoints  0 and  0^: 


x = (x  + x )/2 

0^  n,m  n,m+l 


y0t  = (yn,m  + yn,rH-l)/2 


z = (z  + z . )/2 
n,m  n,nrt-l 


x_  ••  (x  + x )/2 
Oj  n,m  n+l,m 


y0t  = (yn,m  + yn+l  ,tn> /2 


z - (z  + z )/2 

0^  n,m  n+l,m 


(A-4) 
(A-5 ) 
(A-6 ) 
(A-7 ) 
(A— 8 ) 
(A-9) 


This  type  of  receiver  is  found  in  bound  vortex  for  calculating  forces 
as  shown  in  Figure  5. 


Type  3 Receivers 
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The  velocity  is  calculated  at  the  upstream  end  of  a vortex  segment  in 

the  wake  0 : 
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TYPES  OF  SENDERS 

The  vortex  elements  that  induce  velocity  are  called  senders.  Three 
different  types  are  illustrated  here. 


Type  1 Senders 

This  type  of  sender  is  found  in  the  bound  vortex  lattice  for  calcu- 
lating both  circulation  infuence  coefficients  and  forces. 


Coordinates  of  p,q,  and  s are  defined  by 
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Type  2 Senders 

This  type  Is  found  at  the  free  tip  vortex  core 


EDGE  OF 
SIDEWALL 


Type  3 Senders 

This  type  is  found  at  the  free  wake  vortex  lines 
TRAILING  EDGE 

/ S / / S SjT  / / / / * 
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INDUCED  VELOCITY  ON  BOUND  VORTEX  WITH 
UNIT  CIRCULATION 

The  velocity  generated  by  a vortex  segment  is  determined  by  using  the 
Biot-Savart  law.  To  demonstrate  the  calculation,  consider  the  velocity 
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on  a bound  vortex  with  unit  circulation  attributed  from  the  bound  vortex 

itself,  the  tip  vortex  core,  and  the  wake  vortex  lines.  Within  the  bound 

vortex,  the  computation  is  concerned  with  a receiver  of  Type  1.  The  step- 

by-step  procedure  is  as  follows: 

1.  Calculate  the  induced  velocity  due  to  bound  vortex 

a.  The  position  vectors  of  the  control  point  (0  ) of  the 

c 

receiver  with  respect  to  the  ends  of  the  vortex  segment  of  the  sender 
(p»q» s)  are: 


rP  = (xo  - y1  + % • yj  + (z0  ■ yk 

c c r c v 

" (X0  - V1  + (y0  - yq)j  + (z0  ‘ 2q)k 


(A-l 1 ) 


- ye)J  + - zJk 


b.  The  unit  vectors  parallel  to  the  spanwise  and  chordwise  vortex 
segments  of  the  sender  are: 

. _ <».  - V1  * <>„  - y.u  + (2P  - za)i: 

qp  * / o 2 ~2 

V<xp  - v + <yp  - v + <xp  ■ v 


(A-12) 


(xs  - xq)I  + (y,  - ya)3  + (»„  - 

/(xs  - xq)2  + (ys  - yq)2  + (zs  - zq)2 


c.  The  unit  vectors  parallel  to  the  velocities  generated  by  the  span- 
wise  and  chordwise  vortex  segments  are: 


s 

e x r / 

1 e 

vt 

qp  P 

qp 

’ SB 

e x r / 

1 e 

yi 

qs  q 

qs 

(A-13) 


<xp  • V typ ' V <2p " V 

(x„  - x ) (yn  - y ) (zn  - * >i 


x ) + (y  - y ) x (z 
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e x r 
qs  q 


(xb  - V (y8  - V (zs  " Zq) 
(X0  - V (y0  * V (Z0  ‘ Zq> 

^(x  - X)2  + (y  - yn)2  + (zc  - z )2 

8 (|  8 (}  S Cj 


h « | e x r | 
t I qp  pi 


h - ie  x r 

i 1 qs  q 


Here  Equation  (A. 16)  indicates  the  distance  from  the  control  point  to  the 
spanwise  segment  and  Equation  (A. 17)  the  distance  from  the  control  point 
to  the  chordwlse  segment. 

d.  The  unit  normal  vector  of  the  receiver  is  calculated  by  Equations 
(24)  through  (26). 

e.  The  cosines  of  the  angles  between  the  spanwise  segment  and  the 

position  vectors  r and  r are 

P q 

. (A-18) 

(A-19) 

f.  The  cosines  of  the  angles  between  the  chordwlse  segment  and  the 

position  vectors  r and  r are 
q s 


cos  0 


(A-21) 


-e  . r 
qs  a 

‘-s' 

g.  From  the  Biot-Savart  law,  the  velocities  generated  at  the  control 
point  by  the  spanvlse  and  chordwise  segments  (each  have  unit  circulation) 
are 


Vti(D 


(cos  0 + cos  0.)  e 

It 2t  vt 

4wh . 


VU(1)  “ 


(cos  eu  + cos  e2t)  evt 
4irh. 


(A-22) 


h.  The  normal  component  of  the  velocities  at  the  control  point  are 


i 


(A-23) 


i.  Velocity  generated  by  the  image — since  the  upper  end  of  the  side- 
wall  is  attached  to  the  ship,  the  Influence  of  its  flat  upper  surface  is 
accounted  for  by  the  velocity  generated  by  the  image  of  the  sender.  The 
coordinates  of  the  image  are  calculated  by  simply  changing  the  signs  of 
the  z-coordinates  of  the  sender.  For  instance  at  control  point  0c , the 
position  vectors  with  respect  to  p' , q' , s'  of  the  sender  in  the  imaginary 
plane  are 
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V " (x0  " V1  + (yo  ■ yp}j  + (zo  + zP)k 


V ■ (X0  " Xq)I  + (y0  ~ V3  + (zo  + Zq)k 


(A-24) 


V " (X0  ' V1  + (y0  ' yart  + (Z0  + Zs)k 


Following  a similar  procedure,  the  resulting  normal  components  of  the 
velocities  at  control  point  0^  due  to  the  image  are: 


V (2)  - V (2)  . e 
tiij  t j n 


Vtn/2)  " \<2>  * Gn 


(A-25 ) 


j.  The  total  resulting  velocity  generated  at  the  receiver  is  the 
vector  sum  of  those  contributed  by  the  senders  in  the  real  and  image 
planes.  The  normal  components  of  the  velocities  at  the  control  point 
of  the  receiver  (nr)  generated  by  the  sender  (ns)  and  its  image  are  then: 


V (ns.nr)  - V (1)  - V (2) 
tn  till 


(A-26) 


2.  Calculate  the  induced  velocity  due  to  tip  vortex 

a.  The  position  vectors  of  control  points  for  receivers  of  Type  L 
with  respect  to  the  ends  of  tip  free  vortex  are: 


\ ■ (x0  • v1  + (yo  - yj + (z0  ■ zP)k 

v c K c v c v 


rq  " (X0  ‘ Xq)1  + (y0  “ yq)J  + (z0  “ Zq)k 

^ c ^ c ^ c 


(A-27 ) 
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The  unit  vector  parallel  to  the  tip  vortex  segment  is: 


<«„  - »„>*  + <>„  : y°)1  * (z° ~ z°)k 


(A-28 ) 


q ^(x_  - xj2  + (y^  - y„)2  + (zn  - zn)2 


p q 


p q 


p q 


c.  The  unit  vectors  parallel  to  the  velocity  generated  at  the 
control  point  of  Type  1 by  the  tip  vortex  segment  are: 


e x r 

_a 2 E_ 


vt 


(A-29 ) 


e x r I 

qp  P 


h = e x r I 
t qp  p 


(A-30) 


The  second  equation  indicates  the  distance  from  0^  to  the  tip  vortex  seg- 


1 
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d.  The  cosines  of  the  angles  between  the  tip  vortex  segment 


the  position  vectors  r and  r are: 

P q 


cos  e. 


e . r 


(A-31) 


-e  . r 


cos  0, 


(A-32) 


e.  The  velocity  generated  at  the  control  point  by  the  tip  vortt 
segment  with  unit  circulation  is: 


Vti(U 


(cos  8lt  + cos  02t)  ivt 


(A-33) 


f.  The  normal  component  of  the  velocity  at  the  control  point  is: 

Vtn  (1)  " VM(1)  . e 
tnj  tl  n 

g.  The  unit  vector  parallel  to  the  last  tip  vortex  segment  is 


e^p  - (cos  a) i + (sin  a)j 


(A-3A ) 


h.  The  unit  vector  parallel  to  the  velocity  generated  at  control 
point  of  Type  1 by  the  last  tip  vortex  segment  is 


e x r 


I e x r | 

qp  q 


(A-35 ) 


(Here  q is  the  upstream  end  point  of  the  last  tip  vortex.] 
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(A-36) 


h » | e x r | 
t qp  q 

i.  The  cosines  of  the  angles  between  the  last  tip  vortex  segment 
and  the  position  vector  r are 

q 


cos  0 


(x  - x ) cos  a + (y  - y ) sin  a 
0 q "0  q 

c c ^ 


It 


Xq)2  + (y0  " yq)2  + (Z0  " Zq)2 
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cos  0 


2t 
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j.  The  velocity  generated  at  the  control  point  by  the  last  tip 
vortex  segment  is 


V.  (1) 
C1 


(cos  0 + 1)  e 


vt 


4*h 


(A-39 ) 


k.  The  normal  component  of  the  velocity  at  the  control  point 
of  a Type  1 receiver  is 


0 

c 


V (1)  = V (1)  . I (A-40) 

tii  j t^  n 

1.  Velocities  generated  by  the  image  of  the  free  tip  vortex  core 
are  obtained  by  substituting  negative  z-coordinates  for  p and  q in  Equations 
(A-27)  and  (A-28).  The  resulting  normal  component  of  the  velocity  at 
control  point  0^  is 


tn 


(2) 


(2) 


(A-41 ) 
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3.  Calculate  the  induced  velocity  due  to  wake  vortex  segments 

a.  The  position  vectors  of  the  control  point  of  Type  1 receiver 
with  respect  to  a free  wake  vortex  segment  are 


\ ‘ (x0  ' V1  + (yo  ■ V3  + (Z0  " Zq)k 

M c c c 


~TS  ‘ (X0  - XsU  + (y0  ‘ ys)j  + (Z0  ' Zs)k 

c c c 


(A-42) 


q 


b.  The  unit  vector  parallel  to  the  free  wake  vortex  segments  is 


(xs  - x )i  + (ys  - yq)j  - (zs  ~ zq)k 


- *q)2  + (y8  - yq>2  + (*s  - zq>2 


(A-43) 


c.  The  unit  vectors  parallel  to  the  velocity  generated  at  the 
control  point  of  a Type  1 receiver  by  the  free  wake  vortex  segment  are 


e.-e  xr/|e  xr| 
vl  qs  q qs  q 


(A-44) 


h,  ■ |e  x r | 
t 1 qs  q' 


(A-45) 


The  second  equation  indicates  the  distance  from  0c  to  the  wake  vortex 
segment. 

d.  The  cosines  of  the  angles  between  the  free  wake  vortex  seg- 
ments and  the  position  vector  are 


cos  0, 


e . r 
qs  c 


(A-46) 


cos  0, 


-e  . r 
qs  s 


(A-47 ) 


e.  The  velocity  generated  at  the  control  point  of  a Type  1 re- 
ceiver by  the  free  wake  vortex  segment  is 


VU(D 


(cos  0lt  + cos  02l)  evt 


(A-48) 


f.  The  normal  component  of  the  velocity  at  0c  is 


Vln/1)  -\(1)  * en 


(A-49) 


g.  The  induced  velocities  due  to  the  last  wake  vortex  segments 


e - (cos  a)i  + (sin  a)j 
qs 


(A-50) 


e - e x r / 1 r | 
vl  qs  s s 
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(here  "s"  is  the  upstream  end  point  of  the  last  wake  segment.) 
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h,  - | e x r | 
l 1 qs  8 


cos  0,,  - e . r / | r | 
It  qs  8 s 


cos  e2l  - l 


(cos  0U  + cos  e2t)  evt 


V.  (2)  = V.  (2)  . e 
*»n  ^ ^ n 


h.  The  induced  velocities  due  to  the  image  of  the  wake  vortex  are 
obtained  by  substituting  negative  z-coordinates  for  q and  s in  Equations 
(A-42)  to  (A-43).  The  resulting  normal  component  of  the  velocity  at  control 
point  0c  is 


CALCULATION  OF  FORCES  IN  BOUND  VORTEX 

Because  of  length  involvement,  the  velocity  term  in  the  force  equation 
has  to  be  calculated  at  the  middle  point  of  each  vortex  segment  rather  than 
at  the  geometric  center  of  each  panel,  as  in  the  previous  section.  With 
reference  to  the  following  sketch,  therefore,  the  position  vector  of  con- 
trol point  0t  of  a Type  2 receiver  with  respect  to  end  point  p,  of  the 
sender  is,  for  instance: 


51 


I 

1 

I 


!! 


rpt  = (xot  - v1  + (yot  - yj  + (zot  - yk 


Other  position  vectors  r . , r r . , r and  r . assume  similar  forms. 
v pi  qt’  qi  st  si 

The  unit  vectors  parallel  to  the  spanwise  and  chordwise  vortex  seg- 
ments of  the  sender,  e and  e , and  unit  vectors  parallel  to  the  velocity 

qp  qs 

generated  at  midpoints  of  Type  2 senders  by  spanwise  and  chordwise  vortex 

segments  e , e „ , e . , and  e ..  are  then  obtained  by  substituting 
vtt  vlt  vtl  vtl 

Equations  (A-57)  through  (A-62)  for  the  form  of  Equation  (A-13).  Likewise, 
the  velocity  at  a spanwise  receiver  of  Type  2 due  to  a spanwise  sender  is 
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and  the  velocity  at  a spanvlse  receiver  of  Type  2 due  to  a chordwlse 
sender  Is 


(1) 


Yt(ns)  (cos  6Ut  + cos  e2U)  _ 


4"h 


It 


vtt 


(A-59) 


and  the  velocity  at  a chordwlse  receiver  of  Type  2 due  to  a spanvlse  sender 
Is 


(1) 


Yt(ns)  (cos  eut  + cos  e2ti)  _ 

4whtt  evlt 


(A- 60) 


and  the  velocity  at  a chordwlse  receiver  of  Type  2 due  to  a chordwlse  sender 
is 


(1) 


Yj(ns)  (cos  em  + cos  e2tt)  _ 

4whlt  evlt 


(A-61) 


where  the  cosines  are  determined  In  a manner  similar  to  Equations  (A— 18) 
and  (A-19)  with  new  position  vectors. 

Image  contributions  are  obtained  by  changing  the  sign  of  z-coordinates 
and  are  designated  by  "2." 

Finally,  the  velocities  at  the  midpoints  0 and  0,  of  the  receiver 
of  Type  2 by  the  sender  of  Type  1 and  its  image  are 

vot  " vtt(1)  - vtt(2)  + vn(1)  * vu(2)  (A‘62) 

\ " vti(1)  ' vtt(2)  + vti(1)  ' vu(2)  (A_63) 
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The  forces  generated  at  midpoints  0t  and  0^  are  then  calculated  by 
using  the  Kutta-Joukowski  theory 


F - PL  (nr)  y (nr) 

tv  X i ] 

(A-64) 
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t 

F - pL  (nr)  ^(nr) 

lv  x V 

(A-65) 

ut 

t 

where  L (nr)  and  L (nr)  are  the  lengths  of 

spanwise  and 

chordwise  segments 

of  the  receiver,  respectively. 


The  forces  in  bound  vortex  due  to  tip  and  wake  vortices  are  then  cal- 
culated by  following  a similar  procedure. 

EXTENSION  OF  THE  METHOD  TO  DOUBLE- 
SIDEWALL CASES 

The  method  has  been  extended  to  double- sidewall  cases  in  accordance 
with  the  step-by-step  procedure  outlined  below.  Since  the  velocity  and 
force  equations  are  very  similar  in  nature  to  those  presented  previously, 
they  will  not  be  repeated  here. 

1.  Bound  vortex  lattices  for  both  sidewalls  lie  in  x-z  planes. 

The  second  sidewall  is  placed  above  the  first  one  at  y - B. 

2.  The  velocities  with  unknown  circulation  in  bound  vortex  on  Wall  1 
due  to  bound,  tip  and  wake  vortices  of  Walls  1 and  2 are  calculated. 

3.  The  velocities  with  unknown  circulation  in  bound  vortex  on  Wall  2 
due  to  bound,  tip  and  wake  vortices  of  Walls  I and  2 are  calculated. 

4.  The  circulation  values  y are  determined  by  solving  the  influence 
coefficients. 

5.  The  shape  of  the  tip  and  wake  vortex  lines  of  Wall  1 are  then 
updated. 

6.  The  shape  of  tip  and  wake  vortex  lines  of  Wall  2 are  then 
updated. 

7.  Finally,  the  forces  in  bound  vortices  on  Walls  I and  2 are 
calculated. 
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